Many-body quantum dynamics of polarisation squeezing in optical fibre 
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We report new experiments that test quantum dynamical predictions of polarization squeezing 
for ultrashort photonic pulses in a birefringent fibre, including all relevant dissipative effects. This 
exponentially complex many-body problem is solved by means of a stochastic phase-space method. 
The squeezing is calculated and compared to experimental data, resulting in excellent quantitative 
agreement. From the simulations, we identify the physical limits to quantum noise reduction in 
optical fibres. The research represents a significant experimental test of first-principles time-domain 
quantum dynamics in a one-dimensional interacting Bose gas coupled to dissipative reservoirs. 



The nonlinear optical response of standard commu- 
nications fibre provides a straightforward and robust 
method []J for squeezing the quantum noise always 
present in laser light to below the vacuum noise level. 
This feature allows us to design quantum dynamical ex- 
periments 0, 13 operating in a very nonclassical regime 
where highly entangled states can be readily produced, 
even in many-body regimes involving 10 8 interacting par- 
ticles. The squeezing is sensitive to many-body photon- 
photon interactions, as well as to additional dissipative 
and thermal effects |4j . Such complications have affected 
all prior experiments and have so far prevented quanti- 
tative agreement between theory and experiment. 

Here we report on quantitative comparisons of first- 
principles simulations with experimental measurements 
on the propagation of quantum states in optical fiber. 
The excellent agreement, over a wide range of initial con- 
ditions, is unprecedented for direct predictions from ab- 
initio treatments of many-body quantum time-evolution. 
The approach we use has potential applications in many 
other areas of science, especially to dynamical experi- 
ments with ultra-cold atoms and nanotechnology. 

Photons in a nonlinear fibre are an implementation of 
the famous one-dimensional attractive Bose gas model|5j 
in quantum field theory. Fibre squeezing experiments 
thus provide a substantial opportunity to carry out an 
experimental test of the predictions of many-body quan- 
tum mechanics for dynamical time-evolution. Such a 
test requires the same ingredients as did Galileo's famous 
tests of classical dynamics using an inclined plane 0: one 
needs a known initial condition, a well-defined cause of 
dynamical evolution, and accurate measurements. All of 
these essential features are present in our experiments. 
The initial condition is a coherent photonic state pro- 
vided by a well-stabilized pulsed laser. The Kerr non- 
linearity in silica fibre corresponds to a localized (delta- 
function) interaction between the photons Q- Quantum- 
limited phase-sensitive measurements have been devel- 



oped in optics that are able to detect quantum fluctua- 
tions at well below the vacuum noise 

Even though the simplest model of a ID interact- 
ing Bose gas has exactly soluble energy eigenvalues, the 
many-body initial value problem still remains intractably 
complex. Expanding the coherent initial state directly in 
terms of more than 10 100 relevant eigenstates (more than 
the number of particles in the universe) is simply not 
practical. For such reasons, Feynman suggestedfllj that 
the calculation of dynamical quantum time-evolution 
would be impossible without special 'quantum comput- 
ing' hardware, which is not available at present. Even 
worse, there are additional complications occurring due 
to coupling to phonons in the silica fibre medium. 

To deal with the exponential quantum complexity in 
this dissipative system, using existing computers, we em- 
ploy quantum phase-space methods based on the work 
of Wigner and Glauber 0. With modifications to 
treat non-classical fields 0, phase-space methods are 
able to calculate quantum correlations in a dynamical 
system without recourse to hydrodynamic approxima- 
tions or linear response theory. They were successfully 
used for the original prediction of quantum squeezing 
in optical fibres. which was qualitatively verified 
in several laboratories 111 Il5i. Recently, photon- number 
squeezing was simulated and compared with exper- 
iment |17| . However, the cause of squeezing degrada- 
tion was unclear 18j . By using an improved experimen- 
tal configuration and including all relevant dissipative 
effects in the simulations, we have been able to quantita- 
tively test and verify the phase-space predictions. 

Physically, the Kerr effect that generates the squeez- 
ing can be viewed as producing an intensity-dependent 
refractive index, which distorts an initially symmetric 
phase-space distribution into an ellipse Because of energy 
conservation, however, the variance in the amplitude re- 
mains constant. Thus the squeezing cannot be detected 
directly in amplitude or intensity measurements. 
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Figure 1: Experimental detail for the efficient generation of 
polarization squeezed states. Pairs of coherent, 74 fs FWHM 
sech-shaped laser pulses at f.5 fim are coupled into a birefrin- 
gent fiber. A pre-fiber variable delay ensures that the pulses 
overlap upon exitting the fiber. A relative phase shift of ^ is 
maintained by means of a piezo-electric crystal in a feedback 
loop, resulting in circularly or Ss-polarized light. The squeez- 
ing measurement is performed by use of a half-wave plate, po- 
larizing beam splitter and two balanced detectors. Although 
the squeezing is broad-band, it is sampled at a frequency of 
17.5 MHz to avoid technical noise. 



Although laser homodyne measurements can be 
used for these phase-sensitive measurements, they re- 
quire a local oscillator, which is impractical for the 
high-intensity laser pulses used here. Instead, we em- 
ploy an interferometric method known as polarization 
squeezing wherein two multi-photon femptosecond 
pulses that are initially identical in amplitude propagate 
along the mutually orthogonal polarization axes x and 
y of a polarization-maintaining fiber. The two output 
pulses are then combined on a half-wave plate at an an- 
gle $. For the appropriate phase-space rotation angles 
8 = 4$, the squeezing or antisqueezing can be detected 
in the number difference. Because both pulses undergo 
similar evolution, most excess thermal noise is common 
to both modes and is thus cancelled. An exception here 
is any noise induced by birefringent, or depolarizing, ef- 
fects, which will be unique to each mode. The experi- 
mental set-up is illustrated in Fig. 1. 

Because the experiment involves ultrashort pulses, we 
use photon-density operators for each polarisation mode 
of the optical field: $f x (t, z) and ^f y (t, z) that include a 
range of spectral components a(t, k): 



dka a (t,k)e lik " ko)z+luJot , (1) 



where a — x,y labels the polarisation, t is time, z is prop- 
agation distance and k is wave-number. The wavenum- 
ber and frequency of the carrier wave are given by fco and 
luq. The commutation relations of the field operators are 

To describe the polarisation squeezing, we define gen- 



eralized Stokes operators: 

S = N XX (L) + Nyy(L), Si = N XX {L) - N yy {L), 
S 2 = N xy (L) + N yx (L), S 3 = iN yx (L) - iN xy {L)(2) 

where L is the fibre length and N a(T i(z) — 
f dtW a (t, z)\&cr'(t, z). So corresponds to the combined 
intensity, S\ gives the number difference between the x 
and y modes, and S2 and S3 are operators sensitive to 
the relative phase difference. In the Stokes measurement 
after the fiber, we measure a given Stokes parameter or- 
thogonal to the bright excitation in S3. This, the dark 
S1-S2 plane, is scanned by rotating a half-wave plate to 
measure: Se = cos(0)Si + sin(9)S2- Following from the 
Stokes operator commutation relations|l9j. the number 
difference is squeezed if the variance of Se is below, for 
S^-polarized light, 



A 2 S e < 



S'o 



(3) 



To simulate the evolution of the photon density 
^ait, z), we employ a quantum model of a radiation 
field propagating along a fused silica fibre, including the 
^-(3) nonlinear responses of the material and the nonres- 
onant coupling to phonons^, |2£j. The phonons act as 
a reservoir with finite-time correlations (non-Markovian) 
that generates additional, delayed nonlinearity, as well as 
spontaneous and thermal noise. Classical, low-frequency 
phase noise can be neglected, becase it is common to 
both polarization modes, unless there is a depolarising 
component. The phonon spectrum is based on the ex- 
perimentally determined Raman gain for pure fused silica 
Q ;ii ( w )|i3- Because of the fiber birefringence, the two po- 
larization components do not temporally overlap for most 
of the fiber length, so we neglect the cross-polarization 
component of the Raman gain. 

To simplify the theoretical description, we use a prop- 
agative frame with rescaled time r = (t — z/v), rescaled 
propagation distance £ = z/zq and dimensionless photon- 
flux field 4> a = $ a ^Jvts/n. Here to is the pulse duration, 
zq = t^/\k"\ is the dispersion length and 2n is the pho- 
ton number in a fundamental sech(r) soliton pulse. The 
fibre-dependent parameters are pulse velocity v, and dis- 
persion parameter k" . 

From the multimodc Hamiltonian for the full 
system we obtain Heisenberg equations of motion for 
the quantum fields. The equations for the phonon field 
are integrated to derive quantum Langevin equations for 
the photon-flux field 

o - o2 

q^AtX) = ^0a(T,C)+<fa(T,C)&(T,O 



dr'h(T 



r (T',C)^(r',C)? CT (r,C) ,(4) 



where the nonlinear response function h(r) includes both 
the instantaneous electronic response and the Raman re- 
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sponse determined by the gain function q- R fq;')|^.l20l l2l| . 
The correlations of the reservoir fields are 

(ft(o/,c')rv(^c)) = ^^[n th (M) + e(-a,)] 

x«5(C-C')%-^')^'. (5) 

where n t h is the temperature-dependent Bose distribu- 
tion of phonon occupations. The Stokes (ui < 0) and 
anti-Stokes (to > 0) contributions to the Raman noise 
are included by means of the step function O. These 
equations can easily be modified to include coupling to 
gain and absorption reservoirs, but this is unnecessary 
for the short distances in these experiments. 

The exact quantum dynamics can been simulated us- 
ing the positive-P (+P) phase-space representation [Ti) . 
However, for large photon number n and short prop- 
agation distance L, it is known that the +P method 
gives squeezing predictions in agreement with a truncated 
Wigner phase-space method|22j, which allows faster cal- 
culations. We have chosen the latter method to re- 
duce computational overheads. This method maps a 
field operator to a stochastic field: </> ct (C,t) — > (j) a ((, t). 
Stochastic averages involving this field then correspond 
to symmetrically ordered correlations of the quantum 
system. Because of the symmetric-ordering correspon- 
dence, quantum effects enter via vacuum noise. The Kerr 
effect merely amplifies or diminishes this noise in a phase- 
sensitive manner, which makes the Wigner approach ide- 
ally suited for squeezing calculations. 

Using the Wigner mapping, we obtain a Raman- 
modified stochastic nonlinear Schrodinger equation for 
the photon flux that is of exactly the same form as 
Eq. (@}jj,|2fj. The correlations of the Raman noise fields 
To- are also of the same form as their operator equiva- 
lent Eq. except that the step function is replaced 
by a constant 1/2. Because of the symmetrically ordered 
mapping, the Stokes and anti-Stokes contributions to the 
Wigner Raman noise arc identical. The initial condition 
of the stochastic field is the mean coherent level plus fluc- 
tuations that correspond to vacuum noise: 

(A0 ff (r,O)A0;,(r',O)) = -U(r - r')<W • (6) 

In 

Simulations of the pulse propagation reveal the im- 
portance of including the system's multimode nature, 
which affect small-amplitude and intense pulses in differ- 
ent ways. As Fig. 2 shows, the evolution of the amplitude 
profile of a weak pulse is dominated by dispersion. In 
contrast, an intense pulse reshapes into a soliton, whose 
subsequent evolution reveals the effect of the Raman self- 
frequency shift. The range of input pulse energies in the 
experiment includes both these extremes. 

Fig. 3 gives the phase-space rotation angle 9 at which 
squeezing is observed at different input energies. The 
simulations and experiments agree very well for the in- 
tense pulses, but there is a divergence for weak pulses. 




Figure 2: Propagation of (A) E = 4.8 pJ and (B) E = 53.5 pJ 
pulses, with initial width to = 74 fs. 
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Figure 3: Phase-space rotation angle 9 versus total input en- 
ergy E. Squares and diamonds give the experimental results 
for the Li = 13.4 m and L2 = 30 m fibres, respectively. Con- 
tinuous (dashed) lines give the simulation results with (with- 
out) excess phase noise. Other parameters are to = 74 fs, 
zo — 0.52 m, n — 2 x 10 8 , and Ao = 1.51 fim. 

Now classical phase noise has a relatively strong effect 
on the weaker pulses, because the antisqucezing produced 
by the Kerr effect is then smaller and at a greater angle 
to the phase quadrature. We therefore include the effect 
of excess phase noise by a nonlinear least-squares proce- 
dure. A new relative noise variance is calculated from 

A 2 S e /(s 3 J> = p p sm 2 (9 N )+ p s cos 2 (9 N -9 K ) 

+ Pa sin 2 {9 N -6 K ), (?) 

where 9k (E) are the angles for Kerr squeezing with input 
energy E, and where p s (E), p a (E) and p p (E) are the rel- 
ative variances of the Kerr squeezing, Kerr antisqucezing 
and phase noise. The Kerr variances and angles 9k{E) 
are taken from the simulation data, and the phase noise 
is assumed to depend linearly on pulse energy [24]. The 
coefficient of the phase noise is determined by a nonlin- 
ear least-squares fit of the angles 0n(E) that minimise 
the new variance to the measured squeezing angles. The 
result gives an excellent fit to the experimental data, even 
though there is only one parameter for each fibre length. 
This parameter was larger for the longer fiber length, 
indicating that phase noise is largely a fiber-induced ef- 
fect, such as depolarizing guided acoustic wave Brillouin 
scattering (GAWBS). However, there also seems to be a 
smaller bircfringent noise source outside the fiber. 

Antisqueezing and squeezing results are shown Fig. 4, 
for 13.4 m and 30 m of fiber, with the excess phase noise 
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Figure 4: Raw antisqueezing (A) and squeezing (B) for L\ — 
13.4 m (squares) and L2 — 30 m (diamonds) fibers. Error 
bars on the squeezing data indicate the uncertainty in the 
noise measurement; for the antisqueezing the error bars were 
to small to be plotted distinctly for the antisqueezing. Solid 
and dashed lines show the simulation results with additional 
phase noise for L\ and L2, respectively. Dotted lines indicate 
the sampling error in the simulation results. 



included. The theoretical results for both squeezing and 
antisqueezing closely match the experimental data, after 
estimated linear losses of 24% are taken into account. 
That the theoretical antisqueezing results match the ex- 
periment indicates that the fit obtained from the rotation 
angle accounts for all the excess classical noise. A dete- 
rioration of the squeezing is seen at higher intensity due 
to Raman effects, especially for longer fiber lengths. 

As an indication of the accuracy of the simulations, we 
predicted from comparisons of the rotation angle and the 
anti-squeezing curves in theory and in experiment that 
the 30 m fiber should have a 5% larger core diameter 
than the shorter fiber. This was verified by fiber mea- 
surements. There is a residual discrepancy in the highly 
sensitive squeezing measurements for the 30 m case at 
large pulse energies. This could be due to effects such 
as cross-polarization Raman scattering, higher-order dis- 
persion, or initial pulse-shape distortion. 

In conclusion, the efficiency of the squeezing experi- 
ments described here means that a comprehensive theo- 
retical model must be formulated and solved to quanti- 
tatively account for all observations. Solving the model 
requires a first-principles simulation of quantum time- 
evolution in a many-body system coupled to a non- 
Markovian reservoir. We have achieved this by means 
of a Wigner phase-space representation, thus obtaining 
excellent agreement between simulation and experiment 
over a wide range of energies and fibre lengths. The sim- 
ulations reveal that Raman effects limit squeezing at high 
intensity and for longer fibers and that depolarizing phase 
noise (i.e. GAWBS) limits squeezing at low intensity. 
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